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I. INTRODUCTION 

Singly ionized calcium (Ca + ) is an interesting can- 
didate in many areas of physics. It is especially im- 
portant in astrophysics for investigating the radiative 
properties of stellar objects [J 0. Its transition wave- 
lengths and electric dipole amplitudes are required to 
find out isotopic abundances [3J] and the energy trans- 
fers in stars 0, Q . They are also used for obtaining in- 
formation on emission and absorption lines of the elec- 
tric dipole transitions between the low-lying states in 
galaxies, interstellargas clouds and gas disks surround- 
ing the stars [l|, 0, EL @| • Ca + is also suitable for lab- 
oratory physics. Using the techniques of laser cool- 
ing and ion trapping, it has been subjected to many 
precision measurements, optical frequency methodology, 
quantum processing and accurate fine structure constant 
measurements [H 0, H, IH EH EH EH- In these mea- 
surements, the knowledge of polarizabilities is necessary 
to estimate the black-body shift (BBS) and the Stark 
shift due to the external electromagnetic fields. In our 
recent works, we have reported the hypcrfine structure 
constants and quadrupole moments in Ca + using the rel- 
ativistic coupled-cluster (RCC) method [H, E3>EB|- The 
determination of electric dipole polarizabilities requires 
electric dipole (El) matrix elements and excitation en- 
ergies of all the allowed transitions. Due to the impor- 
tance of these quantities, a number of calculations based 
on various many-body methods including the sum-over- 
states approach in the framework of the RCC theory 
are employed to evaluate them [It], El, EH Ull- There 



are also measurements of the static dipole polarizability 
of the ground state in Ca + [^H HH, but the results do 
not agree with each other. In fact, all the calculations 
17, EH EH HH| differ from the recent measured value 
20l | . Therefore, it is necessary to carry out thorough 
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investigations of the role of electron correlation, higher 
order relativistic effects and contributions from the two 
hole-two particle and the neglected one hole- one par- 
ticle excited states in the calculations of polarizabilities 
using an all orderafr initio approach like the RCC theory. 
We have developed a novel technique to account for the 
importance of different correlation effects in these prop- 
erties for closed-shell and one- valence atomic systems by 
directly obtaining the atomic wave functions perturbed 
to first order by the electric dipole operator in the frame- 
work of the RCC theory [1H HH • This method avoids the 
sum-over-states approach and thereby includes different 
types of correlation effects in a rigorous manner. This 
theory has been employed to determine the ground state 
polarizabilities in a few alkali atoms and singly ionized 
alkaline earth ions including Ca + to check the validity 
of the theory [22|]. Although, the theory for the tensor 
polarizabilities has been developed, it has not been ap- 
plied to excited states. In this work, we calculate the 
El matrix elements and excitation energies and employ 
the above approach to determine scalar and tensor po- 
larizabilities of the 4S and 3D states of Ca + . The role 
of the Breit interaction had not been studied in the ear- 
lier works which we investigate here using this ab initio 
method. 

There have been recent measurements on the branch- 
ing ratios from the 4p 2 P^/2 state and the corresponding 
transition probabilities in Ca + [24| which need to be the- 
oretically investigated. We carry out these studies using 
our ab initio approach and by combining our El matrix 



2 



elements with the experimental wavelengths and compare 
with their corresponding experimental results. Wc also 
evaluate the lifetimes of the 4P-states using these results. 

The remaining part of the paper is organized as follows: 
In Sec. II, we present a brief outline of the theory. This 
is followed in Sec. Ill by a discussion of the method to 
evaluate the unperturbed and the first order perturbed 
atomic wave functions using the RCC method. We then 
present the results and discuss the effect of correlation 
on various properties in Sec. IV and in the final section 
we make some concluding remarks. 



II. THEORETICAL APPROACH 

The static dipole polarizability of a state \Jo,Mn) is 
given by 



an 



3M 2 - Jq(J + 1) 
Jq(2J - 1) 



(2.1) 



where and a\ are the scalar and tensor polarizabilities. 
From the angular momentum selection rule, it is obvious 
that «g will be non-zero only for the states with J > 1/2. 
In an explicit form, the expression for the polarizability 
in the sum-over-states approach can be written as 
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and the Es are the energies of the atomic states. In a 
single valence system, a l can be divided into three parts 
in general as follows: 



«o = a l (v) + a l (cv) + a l (c 



(2.3) 



where v, cv and c inside the parenthesis represent for va- 
lence, core- valence and core correlation contributions, re- 
spectively. In the sum-over-states approach, it is custom- 
ary to evaluate a (v) by calculating the important va- 
lence excited states. However, contributions from a Q (cv) 
and a l (c) are generally taken approximately in such an 
approach. On the otherhand, it is possible to calculate 
a exactly in a particular configuration space by evalu- 
ating the wave function that is perturbed by the electric 
dipole operator operator D in the following manner: 



Let us rewrite Eq. 



as 



Eq — E n 



which in Dirac notation can be expressed as 



i nJn - J „ (^ 0) (-/o,7)lig||^ (0 ) (Jn,7 / ))^ (0) (^,7 / )ll^ll^ (0) (Jo,7)) 



En — Er, 



= (*( >(Jo,7)llA||vI/ (1) G/o,7')> + (* (1) (Jo,7')llA||* (0) (^o,7)) (2.5) 

where 7 represents parity eigenvalue of the state | Jo, M ) and 7' is its opposite eigenvalue and we define an effective 
dipole operator as Di = Cj(— 1) Jo ~ J ™D. Here, \^^(Jo, 7')) is the first order perturbation correction to the wave 
function |^°)(J , 7)} due to the dipole operator D and given by 

i*"Vo,y» = Ei^°v...v» ( "'°' (J " 7 :" D ^ m(J - 7 " (2.6) 

It can be equivalently written as 

l* (1) (-/o,7')> = ^-^^I^HJ^y))^ 03 ^^')!^!!^ '^^)) 

±— l* (0) (^^))(* (0) (^^)l^ll* (0) (Jo,7)), 
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(2.7) 
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since the matrix elements between the same parity states vanish. Applying the completeness condition, we get 

(//-S )|* (1) (Jo,7')) - -D\^(J nn )}; (2.8) 



the above equation can be considered as a first order 
perturbation equation arising from D. By solving the 
above equation and Eq. (|2.5p it is possible to evaluate 
Op in the framework of the relativistic coupled-cluster 
theory. 
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FIG. 1: Schematic low-lying energy level diagrams and decay 
channels of the P-states in Ca + . 

The branching ratio (BR) of a state (/) to a lower 
energy state (i) is defined as 



-4, 



(2.9) 



where Af^i is the transition probability of the corre- 
sponding transition and sum over i represents total prob- 
abilities of all possible transitions. As shown the low- 
lying energy levels of Ca + in Fig. [TJ the electrons from 
the Ap 2 P\/2 state will jump either to the 4s 2 5 , 1 / 2 or 



3d 2 D 3 / 2 states due to the allowed transition with dif- 
ferent probabilities. Again, electrons from the 4p 2 P 3 /2 
state will jump to the 4s 2 S , 1 / 2 , 3d 2 D 3 / 2 an d 3d 2 D 5 / 2 
states due to the allowed transitions and 4p 2 P\/2 state 
due to the Ml forbidden transition. The lifetime of the 
state can be determined from 

Tf = ■ (2-10) 



By combining Eq. 

A 



and Eq. ([2~TU|) . it yields 



f- 



T f r f- 



(2.11) 



From the same or different measurements of Tf and T f^i, 
it is possible to estimate the corresponding A f—>i for var- 
ious transitions. 

The probabilities due to El and Ml transitions are 
given by 



A E1 



2.02613 x 10 18 



(2J/ 
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and 



A M1 



2.69735 x 10 13 
(2J> + 1)A 3 ^ 



\(J f \\Mm)\ 2 , (2.13) 



respectively. In the above equations, \f—,i is the wave- 
length of the corresponding transition and it is the recip- 
rocal of the excitation energy (EE). 



III. METHOD OF CALCULATIONS 

The RCC method which is equivalent to all order per- 
turbation theory has been recently used to obtain precise 
results and account for the correlation effects in single 
valence systems [HI, 0, [HI . Atomic wave functions for 
single valence systems can be expressed in the framework 
of RCC theory as 



l^ 0) ) = =e T {l + S v }\^ 



(3.1) 



where |$„) is the reference state constructed from the 
Dirac-Fock wave function |$o) of the closed-shell config- 
uration ls 2 2s 2 2p 6 3s 2 3j? 6 }(= |[3p 6 ])) of Ca+ by defining 
= aJ|$o) with a J representing addition of a valence 
electron v. Here T and S v are the RCC excitation oper- 
ators which excite electrons from | [3p 6 ] } and aj|[3p 6 ]) for 
the corresponding valence electron v, respectively. The 
amplitudes of these excitation operators are solved by 



{$ L \{H N e T }\* ) =0 



(^\{H N e T }S v \$ v ) 



-{*i\{H N eT}\* v ) 
+{^\S V \9 V )AE V , 



(3.2) 



(3.3) 



with the superscript L(— 1, 2) representing the single and 
double excited states from the corresponding reference 
states and the wide-hat symbol over H^e T represent the 
linked terms of normal order atomic Hamiltonian iJ/v and 
RCC operator T. AE V is the corresponding valence elec- 
tron affinity (negative of the ionization potential (IP)) 
energy which is evaluated by 



AE V = (<S> v \{H N e T }{l + S v }\$ v ) 



(3.4) 



/- 



The EE between two different states are determined from 
the difference of their AE v s. In Eqs. (|3.2[) and (|3.3p we 
have considered only the single and double excitations, 
however we have incorporated contributions from impor- 
tant triple excitations to the AE V calculations. After 
obtaining the amplitudes for T, the core excitation op- 
erator, we solve Eqs. (|3.3p and (|3.4p simultaneously to 
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TABLE I: Ionization potentials (in au) of Ca from different 
works. 



State 


This work 


Others 


Expt. [25] 


4s 




-0.43627757 


-0.43836° 


-0.43627767 






-0.43802 ft 
-0.436287 c 




3d 




-0.37396663 


-0.37407 a 


-0.37408278 






-0.37485 6 
-0.373921 c 




3d 


2 D 5/2 


-0.37361011 


-0.37379 


-0.37380626 






-0.37448 6 
-0.373921 c 




4p 


2 t-> 

-M./2 


-0.32123958 


-0.32217 a 


-0.32149667 






— 0.32224 6 
-0.320844 c 




4p 




-0.32025203 


-0.32111" 


-0.32048108 






-0.32118 6 
-0.320844 c 




5s 


2 S\/2 


-0.19788645 


-0.198293 c 


-0.19858760 


4d 


D 3/2 


-0.17674718 


-0.175144 c 


-0.17729894 


4d 


2 D i/2 


-0.17665912 


-0.175144 c 


-0.17721142 


5p 


2 Pl/2 


-0.15978307 


-0.160060 c 


-0.16046888 


5p 


2 P 3 /2 


-0.15944054 


-0.160060 c 


-0.16011231 



"Relativistic MBPT(2)[2e_ 
^Brueckner approximation [22^ 
c Non-relativistic Coulomb approximation [l{ 



obtain the amplitudes for the S v operator. We use the 
Dirac-Coulomb-Breit Hamiltonian which is given by 



state in the presence of the electric dipole operator D as 

|C) = = e T+n {l + S v +A v }\$ v ), (3.6) 

where f2 and A v are the modified RCC operators to the 
T and S v operators, respectively. Since Eq. (|2.8p is first 
order in the D operator, the above expression will reduce 
to 

|*„) - = e T {l + S v + fl(l + S v )+A v }\<S> v ). (3.7) 

Now, separating the above wave function as \^ v °^) and 
we get 

- =e T {n(l + S v ) + A v }\$ v ). (3.8) 
Following Eq. (|2.8[) , we solve again the amplitudes for 
the modified operators as 

($ L |{5^fi}|0> ) = -(<t> L \De^\%) (3.9) 

(1 + S v )}\$ v ) + ($5\A V \$ V )AE V , 

(3.10) 

where De T represent again the connecting terms between 
D and T operators. In the singles and doubles approxi- 
mation, we write 



T 

n 



T!+T 2 

fix + n 2 



(3.11) 
(3.12) 
(3.13) 



H = ca • p + {13 - l)c 2 + V nuc {r) + 
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(3.5) 



where c is the velocity of light, a and (3 are the Dirac 
matrices and V nuc {r) is the nuclear potential. 

We extend the RCC ansatz for the perturbed atomic 



and 



A, 



Aiu + A 2l 



(3.14) 



where the subscripts {1,2} represent the single and dou- 
ble excitations, respectively. 

Now the expression for the dipole polarizability follows 

as 



J 



<¥ v 0) \¥ v 0] > 

+ + S v ) + A„}|$„) + ($„|{At + (1 + 5t)nt}pT{i + s„}|$„) 

{1 + S V }N {1 + S v } 



(3.15) 



where we define L>i — (e T DiC T ) and No = e T e T . Gen- 
erally, both Di and iVo in the RCC approach are each 
represented by a non-terminating series. However, we 



have devised a procedure motivated by physical consid- 
erations to deal with them using the Wick's generalized 
theorem. We evaluate first the effective zero-body, one- 
body, two body terms etc. systematically and then sand- 
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TABLE II: Transition matrix elements (in au) from different 
calculations. Recommended values from our work is given as 
reco. 



Transition 


STOs 


This work 
GTOs 


reco 


Others 


4p 2 Pi/ 2 -> 4s 2 S 1/2 


2.86 


2.90 


2.88 


2.890° 










2. see 6 -' 










2.861°'" 










2.898 c 


4p 2 Pi/ 2 -> 3d 2 D 3/2 


2.50 


2.41 


2.40 


2.373° 








2.410°'' 










2.244°'" 


4p 2 -P 3 / 2 -> 4s 2 Si/ 2 


4.02 


4.09 


4.03 


4.088° 








4.060°'' 










4.059°'" 










4.099 c 


4p 2 P 3 /2 -» 4p 2 P 1/2 


1.15 


1.15 


1.15 




4p 2 P 3 /2 -> 3d 2 D 3/2 


1.12 


1.09 


1.09 


1.059° 








1.076°'' 










1.028°'" 


4p 2 P 3 /2 -> 3d 2 D 5/2 


3.36 


3.28 


3.22 


3.186° 








3.234°'' 










2.995°'" 










3.306 c 


"Relativistic MBPT(2)|2f| 
6,1 Length gauge result with Brueckner approximation 12711 
''"Velocity gauge result with Brueckner approximation [27j| 
c Linearized RCC method [l7t| . 

wich them (except zero-body terms) between the S v , A v 



and their conjugate operators. We have successfully ap- 
plied this method in our earlier works [13, LLJ, [la, [22, [23| • 
The above zero-body terms, open-terms connecting only 
with fl and terms with A„ give us core (ag(c)), core- 
valence (c^q(cv)) and valence (a (vj) correlation effects, 
respectively. 

We also explicitly present contributions from the nor- 
malization factors evaluating them in the following way 



where i represents the grid points which we have taken as 
750 in total, the step size h is taken as 0.03 in the present 
case and ro is the starting point of the radial distribution 
from where the electron orbitals become finite and taken 
as 2 x 10~ 6 . The STOs and GTOs are given by 



Norm = 



l 



l + N v 



-1}, 
(3.16) 



where JV, = {1 + S$}JV {1 + S v }. 



IV. RESULTS AND DISCUSSIONS 

We have employed two different types of the basis func- 
tions to generate the atomic orbitals; Slater type orbitals 
(STOs) and Gaussian type orbitals (GTOs). These or- 
bitals arc defined on a grid given by 



oKi-i) 



1 



F 



STO 



in) 



and 



7 GTO 



(n) 



(4.2) 



(4.3) 



respectively. Here n K is the radial quantum number of 
the orbitals and on is a parameter whose value is chosen 
to obtain orbitals with proper behavior inside and outside 
the nucleus of an atomic system. We further define on as 



(4.4) 



(4.1) 



We have considered a = 0.0975 and (3 = 1.77 for STOs 
and a = 0.00525 and = 2.83 for GTOs. However, 
we have taken 35, 35, 30, 30 and 25 STO and GTO ba- 
sis functions to construct the s, p, d, f and g orbitals 
respectively. For RCC calculations, we have considered 
all the core orbitals and virtual orbitals are considered 
up to 3500 au for s, p and d symmetries and 1500 au 
for f and g symmetries in the present calculations. In 
fact, it is observed that number of virtual orbitals ob- 
tained using STOs are more in a given upper energy limit 
than GTOs while bound orbital energies match well in 
both the cases. To account for the contributions from the 
high lying orbitals in some of the properties that we have 
considered, we have estimated contributions from virtual 
orbitals using the second order many-body perturbation 
theory (MBPT(2)) and recommended (reco) results are 
given by taking into account all these contributions. 

In Table [J we present our IP results for the low- lying 
states and compare them with the corresponding exper- 
imental results. These results using STOs and GTOs 
were consistent. Some IPs from the excited states devi- 
ate from the experimental results and it might be pos- 
sible to improve them by increasing the virtual space. 
We also compare our results with other theoretical re- 
sults. Guet and Johnson had employed the relativistic 
MBPT(2) method to obtain their results [2f|. Liaw had 
employed the Brueckner approximation method to eval- 
uate these energies [27j and his results match with the 
above MBPT(2) results. In a recent work, Mitroy and 
Zhang have used a one electron semi-empirical core po- 
tential in the non-relativistic framework [19|] to estimate 
these energies which cannot distinguish the fine structure 
levels. Our method in contrast is ab initio and electron 
correlation effects are included to all orders in perturba- 
tion theory in the residual Coulomb and Breit interaction 
in the one hole-one particle, two hole- two particle and 
partial three hole-three particle approximation. 

We present the El and Ml matrix elements in Table HT1 
As can be seen, results from our STOs and GTOs differ 
for different transitions. We have considered contribu- 
tions from virtual orbitals from both the basis functions 
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TABLE III 


: Transition probabilities 


(in xlO 6 


s- 1 ) in Ca+. 


Transition 




This work 


Others 












4P *i/2 — * 


4s 2 S 1/2 


135.240 


134.333 


lOO. ZD 






132. 9 b 




















loO.U 


Ar. 2 O 

4p Pl/2 — * 


3d 2 D 3/2 


9.0431 


8.971 


0. / / 






9.0 










1.0 










Q /I 


Ar* ^ p 

4p -T3/2 — * 


4s 2 S 1/2 


135.842 


135.036 


loo.yo 






136. 9 b 










loD.y 










1 on 7^ 
139. f 


4p 2 P 3 /2 -» 


4p 2 Pi/ 2 


~ io- 10 


~ IO" 10 




4p -H3/2 


3d 2 D a/2 


1.055 


0.962 


u.yo 










0.95 










0.87 c 










0.997^ 


4p 2 P 3 /2 -» 


3d 2 D 5/2 


8.435 


8.419 


8.24 a 






8.5 b 










7.2 C 










8.877 d 



"Relativistic MBPT(2) is used [20] 
'Length gauge result with Brueckner approximation [27} 
c Velocity gauge result with Brueckner approximation [27 
d Linearized RCC method is employed (l7|. 



using MBPT(2) and finally given the consistent results 
as reco values. Guet and Johnson [26] have used B-spline 
basis based MBPT to obtain these results. Again, Arora 
et al. [17| have also used a B-spline basis but a linearized 
RCC method to obtain their results . Liaw (2?| has used 
the Brueckner approximation method to get El matrix 
elements in both the length and velocity gauge expres- 
sions. Our method intrinsically contains all these many- 
body effects. We have also evaluated Ml matrix element 



between the 4p 
around 1.15 au: 



'P 



3/2 



4p P±/2 transition which is 



almost same with the 3d D 



5/2 



3d 



2 D 3 / 2 transition [291 ]. 

Using the above matrix elements, we determine the 
transition probabilities and present them in Table IIIII 
We have followed two approaches to determine them. 
First we have considered energies from our calculations 
and derived wavelengths (A ca ') to obtain the ab initio re- 
sults. In the other case, we use our matrix elements with 
the experimental wavelengths (X expt ). Although the Ml 
transition amplitude from the 4p 2 P 3 /2 state is finite due 
to a very small fine structure splitting, the correspond- 
ing transition probability is almost negligible. We have 
also compared our results with other ab initio and semi- 
empirical results in the same table. In a recent work, 



Gerritsma et al. [2J] have measured BRs (we discuss 
these results below in detail) from the 4p 2 P%/2 state and 
obtain various transition probabilities from this state by 
combining their results with the lifetime measurements 
as given by Eq. (|2. 1 If) . Our results with \ ex v l match well 
with their results. 

Using the above transition probabilities, we determine 
BRs from different calculations and present them in Ta- 
ble IIVI These results are compared with the recently 
measured values of the 4p 2 P 3 / 2 state [24j |, As presented 
in this table, our results with A exp * match well the mea- 
surements. When we evaluate BRs for the 4p 2 P\/2 and 
4p 2 P^/2 states due to 3d states using the relation 



i=3d 2 D 3 /2:3d 2fl 5/2 A f->i 



(4.5) 



it gives as 14.97 and 14.40, respectively, which are not 
within the error bar of the existing experimental result 
[28| and hence require further measurements for verifica- 
tion. 

There are a number of experimental lifetime measure- 
ments available for the 4p 2 Pi/2 and 4p 2 P^/2 states 

using beam laser, beam foil, 
beam foil with cascade correction and Hanle techniques. 
Among them the laser-beam-ion-beam spectroscopy by 
Jin and Church [3(| results are the most precise. Sub- 
stituting our transition probabilities in Eq. (|2.10j) . we 
obtain the lifetimes of the 4p 2 Pi/ 2 and 4p 2 P 3 / 2 states 
as 6.931s and 6.881s with A ca ', respectively, where as 
6.979s and 6.924s with A ea;pt , respectively. Other calcu- 
lations based on the above discussed results also predict 
results close to ours. In fact, our result 6.924s of life- 
time of the 4p 2 P-ij2 state is in good agreement with the 
experimental results. 

Using the same wave functions used to obtain the 
above properties and solving Eq. (|2.8|) . we obtain the 
static dipole polarizabilities of the 4s 2 5i/2, 3d 2 D 3 /2 
and 3d 2 D 5 / 2 states with STOs and GTOs and they are 
presented in Table fVTl The dipole polarizabilities for the 
ground state from STOs and GTOs are in good agree- 
ment, but the 3d state dipole polarizabilities differ by 
4%. Since we were able to generate less number of vir- 
tuals using GTOs in a given energy upper bound, the 
convergence of these results were checked with virtual 
orbitals with higher energies which was not possible for 
STOs due to the computational limitation. Therefore, 
we consider our results from GTOs are more accurate 
than results from STOs. There are also a number of the- 
oretical calculations available on both the ground and 3d 
excited states inclu ding our previous work and references 
therein [U OJ, [H H [37j]. We had just carried out 
the ground state polarizability calculation in Ca + along 
with other atomic systems in the earlier work [22j to ver- 
ify the validity of the method that was proposed for the 
first time. In the present case, we have investigated the 
accuracy of the wave functions in Ca + to obtain IPs and 
El matrix elements which are the ingredients to evaluate 
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TABLE IV: BRs of 4p 2 P 1/2 and 4p 2 P 3/2 states in Ca+. 



Transition 






This work 


Others 


Expt [24] 












a 2 n 

4p 2 A/2 -> 


a 2 n 

4s 2 S 1/2 


0.9373 


0.9374 


0.9391° 








0.9366° 
0.9444 c 
0.9350 d 




4p 2 Pi /2 -» 


3d 2 D 3 / 2 


0.0627 


0.0626 


0.0609° 








0.0634° 
0.0556 c 
0.0650 d 




4p 2 P 3/2 -» 


4s 2 S 1/2 


0.9347 


0.9350 


0.9381° 
0.9354° 
0.9443 c 
0.9340 d 
0.9357 6 


0.9347(3) 


4p 2 P 3/2 -> 


4p 2 A/2 


~ 








4p 2 P 3 /2 -> 


3d 2 D 3/2 


0.00726 


0.00666 


0.00628° 


0.00661(4) 






0.00649° 
0.00600 c 




4p 2 P 3 /2 


3d 2 D 5/2 


0.0581 


0.0583 


0.0556° 
0.0581° 


0.0587(2) 



0.0497 c 
0.0593 d 
0.0643 e 



°Relativistic MBPT(2) is used [20] 
"Length gauge result with MCDF method [23 
c Velocity gauge result with MCDF meth od 127 
^Linearized RCC method is employed [l7| ] 
e Semi empirical pj| 



accurate dipole polarizabilitics. In fact, the correlation 
behavior in the 3d state dipole polarizabilities is not dis- 
cussed in the literature. Patil and Tang [36] had used 
multipolar-matrix elements based in the non-relativistic 
approximation to obtain the 4s 2 5 , 1 / 2 state dipole polariz- 
ability. This has got both the summation and integration 
approach over the intermediate states from different or- 
bital quantum numbers. Using Coulomb approximation 
with the Hartree-Slater core calculations, Theodosiou et 
al. [13] had reported the dipole polarizability of the 
same state. Their result differs from ours and it seems 
as though they have not taken core-correlation into ac- 
count. Recently, Arora et al. [13] and Mitroy and Zhang 
[lj| have also evaluated dipole polarizabilities based on 
the sum-over El matrix elements and oscillator strengths 
between different states. The main differences in their re- 
sults and ours is that they have estimated core (neglected 
for tensor polarizability) and core-valence correlation ef- 
fects approximately whereas we have used the first order 
perturbed RCC method to evaluate them. Contributions 
from the continuum and doubly excited states with con- 
figurations like [Ap b ]nsms [n ^ m, with n, m being prin- 



cipal quantum numbers) which are also important for the 
dipole polarizability calculations of the states have been 
considered by us. They are implicitly accounted for in 
the present work by evaluating the first order perturbed 
wave functions due to the electric dipole operator. We 
have also corrected our results due to the normalization 
of the wave functions. In Table IVII1 we present contri- 
butions from the DF and the individual RCC terms ob- 
tained using GTOs. The differences between these two 
results give the correlation contributions associated in 
evaluating these quantities. It is evident from our stud- 
ies that correlation effects in the 3d-states are more than 
50% while it is about 20% in the 4s 2 S 1/2 state. The 
olq{c) and a^cv) contributions are found to be smaller 
for the scalar dipole polarizability than the previously 
estimated results. We also present these contributions 
for the tensor polarizabilities which were neglected ear- 
lier. Contributions due to the doubly excited states and 
normalization corrections cannot be neglected in preci- 
sion calculations. There are three experimental results 
for the ground state dipole polarizability [10, HH , but 
they do do not match with each other. Although the re- 
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TABLE V: Lifetimes (in s) of 4p 2 P 1/2 and 4p 2 P 3 / 2 states in 
Ca+. 



and Zhang 1S| 
0.39(27)Hz 



State 


This work 


Others 


Expt 


4p 2 Pi/ 2 


6.931 6.978 


6.94 a 


7.098(20) d 






7.047 M 


7.07(7) e 






7.128 M 


7.5(5/ 



4p 2 P 3 /2 6.881 6.924 



6.875 c 
6.75 a 
6.833 M 

6.898 M 

6.686 c 



6.62(35) 9 

6.924(19) d 
6.87(6) e 
7.4(6) / 

6.68(35) 9 
6.72(2)'' 

6.61(30)* 



"Relativists MBPT(2)[2J| 
b ' 1 Length gauge result with MCDF method [27|l 
^Velocity gauge result with MCDF method [2^] 
c Linearized RCC method [13]. d Laser-beam-ion-beam 
technique 30] 
e Laser-beam technique s 13111 
' Beam foil technique [32] 
9 Beam foil technique with cascade correction [53] 
h Hanle method [341 
1 Hanle method 351 



suit given by Theodosiou et al. [20( is the latest, but our 
results are close to Chang [2l| . 

The frequency shift (in Hz) due to a black-body (BBS) 
due to the frequency-dependent electric field at temper- 
ature T = 300K by neglecting the dynamic correction in 
the 4s 2 S*i/2 — * 3d 2 D 5 /2 transition is approximated by 

Si 



Au = i(831.9V/m) 2 [a J (4s) _ a J(3d 5/2 )f4.6) 

By substituting our results in the above expression, 
we obtain Av = 0.376Hz which is in agreement with 
0.38(l)Hz by Arora et al. 17] and 0.368Hz by Mitroy 



This also supports the measured value 
10| . The agreement between different calcu- 
lations is mainly due to the cancellation of the results of 
4s 2 Si/ 2 and 3d 2 D 5 / 2 states. 

In Table IVIII1 we present the contributions from the 
Breit interaction to different properties. These contribu- 
tions are smaller in these properties than in the hyperfine 
structure constants which were reported recently [14| . In 
contrast to the hyperhne constants where the Breit inter- 
action contributes more to the 4s 2 Si/ 2 state, it is larger 
in the 3d states than the ground state in the these prop- 
erties. 

V. CONCLUSION 

We have employed the relativistic coupled-cluster 
method with two different basis functions to study ion- 
ization potentials, electric dipole matrix elements and 
dipole polarizabilities in the singly ionized calcium. We 
have also evaluated transition probabilities, branching ra- 
tios and lifetimes of the first excited p-states using these 
results. By determining the first order perturbed wave 
function due to the electric dipole operator, we obtain ab 
initio results for the static dipole polarizabilities in the 
ground and first excited d-states. Black-body shift in the 
4s 2 Si/ 2 — > 3d 2 D 5 / 2 transition has been evaluated us- 
ing these results and compared with the other available 
results. Contributions from the Breit interaction to the 
above properties have been studied for the first time in 
singly ionized calcium. 
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State 




This work 






Others 


Expt 




GTOs 


(reco) STOs 












ah 






a 






73.002 


74.342 




(0.1(1.1 ) 




7H on/1 c\e 

r U.oy^lo J 










(0.49 




Q f A\ f 
fQ.3{4) J 










7fl 070C 

(0.0(2 




70 r / 1 ^^ <? 










70.6 d 






3d 2 D 3/2 
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-15.870 31.604 
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32.73 c 


-25.20 c 










25.4 d 
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32.0(1.1) 
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-25.20 c 
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Terms 
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3d 2 D 3/2 




3d 2 D 5/2 


DF 


96.201 
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-0.268 
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Si v DiAi v + cc 
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-5.161 


4.277 


S2vDiA 2v + cc 
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-0.298 


0.105 


-0.289 


0.132 


Others 


0.130 


0.161 


-0.226 


0.161 


-0.050 


Norm 


-0.754 
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Note: Subscripts lv and 2v represent the valence contributions due to the singly and doubly excited states, respectively. 
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